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Abstract 

The magnetic order-disorder layering transitions of a spin-1 Ising model are investi- 
gated, under the effect of a variable surface crystal field A^, using the mean field theory. 
Each layer k, of the film formed with N layers, disorders at a finite surface crystal field 
distributed according to the law Ak = A^/fc", k = 1,2,...,N and a being a positive 
constant. We have established the temperature-crystal field phase diagrams and found a 
constant tricritical point and a reentrant phenomenon for the first /cq layers. This reen- 
trant phenomenon is absent for the remaining N — k^ layers, but the tricritical points 
subsist and depend not only on the film thickness but also on the exponent a. On the 
other hand, the thermal behaviour of the surface magnetisation for a fixed value of the 
surface crystal field A^ and selected values of the parameter a are established. 
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1 Introduction 



The Blume-Capel model (BC) was originally proposed to study the first-order magnetic phase 
transitions in spin— 1 Ising systems [1]. This model was generalised to the Blume- Emery- 
Griffiths (BEG) to study phase separation and superfluity in ^He-^He mixtures [2]. Later it 
has been applied to describe properties of multicomponent fluids [3], semiconductor alloys [4] 
and electronic conduction models [5]. The (BC) model is not exactly solvable in more than 
one dimension, but it has been studied over infinite (/-dimensional lattices by means of many 
different approximate techniques and its phase diagram is well known. 

When the theory of surface critical phenomena started developing, some attention has been de- 
voted to the study of the (BC) model over semi-infinite lattices, with modified surface couphngs. 
Benyoussef et al. [6] have determined the phase diagram in the mean field approximation, re- 
porting four possible topologies at fixed bulk/surface coupling ratios. A similar analysis have 
also been done using a real space renormalization group transformation [7]. Other works re- 
ferring to particular regions of the phase space are those using: the mean field approximation 
[8], the effective field approximation [9] and the low temperature expansion [10]. All these 
works show that it is possible to have a phase with ordered surface and disordered bulk, which 
is separated from the completely ordered phase by the so-called extraordinary transition and 
from the completely disordered phase by the surface transition. When such a phase is absent, 
the transition between the completely ordered and the completely disordered phase is called 
ordinary. The meeting point of the lines of these three kinds of phase transitions is named 
special and it is generally a multicrititcal point. As discussed in Ref. [11], the strong interest in 
these models arises partly from the unusually rich phase transition behaviour they display as 
their interaction parameters are varied, and partly from their many possible applications. In 
most of these cases considered so far the bilinear interaction is ferromagnetic. In the antiferro- 
magnetic case, the spin-1 Ising systems are used to describe both the order-disorder transition 
and the crystallisation of the binary allow, and it was solved in the mean field approach [12]. 
One of the most interesting and elusive features of the mean field phase diagram for the anti- 
ferromagnetic spin-1 Blume-Capel model in an external magnetic field is the decomposition of 
a line of tricritical points into a line of critical end points and one of double critical points [13]. 
This model was also studied by transfer-matrix and Monte Carlo finite-size-scahng methods 
[14] , but such decomposition does not occur in this two dimensional model. On the other hand. 
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ferroelectric films can be described by an Ising model and when the film becomes very thick, 
its properties are those of the semi-infinite Ising system [15-17]. From the experimental point 
of view, the most commonly studied magnetic multilayers are those of ferromagnetic transi- 
tion metal such as Fe/Ni, where the couphng can exist between magnetic layers [18-20]. The 
discovery of enormous values of magnetoresistance in magnetic multilayers are far exceeding 
those found in single layer films and so exceeds the discovery of oscillatory interlayer coupling 
in transition metal multilayers. These experimental studies have motivated much theoretical 
works to study magnetic thin films as well as critical phenomena [21-26]. This is partly moti- 
vated by the development of new growth and characterisation techniques, but perhaps more so 
by the discovery of many exciting new properties, some quite unanticipated. Using the mean 
field theory, Benyoussef et al. [27] and Boccara et al. [28] have studied the spin-1 Ising model 
with a random crystal field. 

The effect of the surface and bulk transverse fields on the phase diagrams of a semi infinite 
spin-1 ferromagnetic Ising model with a crystal field was investigated in [29] within a finite 
cluster approximation with an expansion technique for cluster identities of spin-1 localised spin 
systems. On the other hand, the transverse field or crystal field effects of spin-1 Ising model 
has been studied by several authors [30-33]. 

The purpose of this paper is to study the effect of a variable crystal field according to the law 
Afe = As/A;° (As being the surface crystal field and k the layer number from the surface and a 
a positive constant), on the order-disorder layering transitions of a spin-1 using the mean field 
theory. This paper is organised as follows. Section 2 describes the model and the method. In 
section 3 we present results and discussions. 

2 Model and method 

The experimental measurements of layer-by- layer ordering phenomena have been established on 
free-standing liquid crystals films such as nmOBC (n-alkyl-4'-n-alkyloxybiphenyl-4-carboxylate) 
[34,35] and bACOOBC (n-pentyl-4'-n-pentanoyloxy-biphenyl-4-carboxylate) [36] for several 
molecular layers. More recently, Lin el al. [37] have used the three-level Potts model to show the 
existence of layer-by- layer ordering of ultra thin liquid crystal films of free-standing 5AC00BC 
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films, by adjusting tlic inter layer and intralayer couplings between nearest-neighbouring molecules. 
The system we are studying here is formed with N coupled ferromagnetic square layers in the 
presence of a crystal field. The Hamiltonian governing this system is given by 

E J^JS^SJ+Y^\{S,f (1) 

<i,j> i 

where, Si{l = = —1,0, +1 arc the spin variables. The interactions between different spins 
are assumed to be constant so that Jij — J. The crystal field acting on a site i is so that 
Aj = Afe for all spins of the layer k so that: Ai > A2 > ... > Ajv- 

Using the mean field theory, the quadrupolar moment and the magnetisation of a plane k are 
given, respectively, by : 

qk -< [bk) >- 1 + 2 cosh(/3/ifc) exp(-/3Afc) ^ ^ 

2sinh(/3/ifc)exp(-/?Afc) 
1 + 2 cosh{fjhk) exp(— pAjt) 

where hk — J{^mk + ruk+i + and (5 — l/{kBT), ks the Boltzmann constant ant T the 

temperature. 

The model we are studying, in this paper, corresponds to a crystal field distributed according 
to the law: 

A, = As/k- (4) 

where A^ = Ai is to the crystal field acting on the surface (first layer k — 1). Values of the 
parameter a will be discussed in the following; in particular a = corresponds to a uniform 
crystal field A^ applied on each layer of the film. 
The free energy of a layer /c, can be expressed as: 

Fk = log(l + 2 cosh(/3/ifc)) + ^mk{imk + m^+i + m^.i) (5) 

the parameters hk and P still have the definitions given above. 

As far as we know there is no natural compounds with decreasing crystal field from the 
surface to deeper layers. However, there exist a series of single crystals R2FeuB {R — 
Y, Nd, Gd, Th, Er, Tm) with a variable crystal field anisotropy. Indeed, by adjusting this 
anisotropy for each layer so that the amplitude of the crystal field will decrease from a layer k 
to the next layer + 1. Hence, these materials can be grown layer-by- layer in a thin film with 
a decreasing crystal field from the surface to deeper layers. This can lead to an experimental 



realisation of our model. Moreover, the decreasing function considered in this model, i.e. k~°', 
is not unique since the order-disorder layering transitions are found for any decreasing function 
of the crystal field from the surface to deeper layers. 

3 Results and discussion 

We are considering a system formed with ferromagnetic square layers of a simple cubic 
lattice spin— 1 Ising model, with free boundary conditions. The notation D'^O^'^ will be used 
to denote that the first k layers are disordered while the remaining N — k layers are ordered. 
In particular, corresponds to an ordered film whereas denotes a totally disordered film. 
The surface crystal field A^, applied on each layer k is distributed according to the law given 
by Eq. (4). 

The ground state phase diagram of this system is illustrated by Fig. 1. For very small values 
of the surface crystal field A^, the system orders in the phase O^. When increasing A^ the 
surface (first layer k — 1) disorders and the phase DO^~^ occurs at A^/J = 3(1)"^. Increasing 
As more and more, the second layer k = 2 becomes disordered at A^/J = 3(2)", and so on. 
The transition from the phase D'^Q^^'^ to the phase is seen at A^/J = 3{k + 1)" 

provided that k + 1 < N. For higher values of the surface crystal field the system is totally 
disordered and the phase occurs. 

In particular, it is found that there exist a critical order layer ko corresponding to the transition 
]jkoQN-ko ^ jjN^ ^ reduced surface crystal field given by: 

N 

A,/ J = (3(7V - A;o) - !)/( ^ (1/n)- (6) 

fe=feo+l 

ko is exactly the number of layering transitions existing at T = 0. It is shown that ko depends 
both on the parameter a and the film thickness N. The special case: o; = is a situation with 
a constant crystal field apphed on each layer and there is only a single transition D^, 
occurring at Ag/ J = 3 — 1/N for T = 0. 

A preliminary study, for N — 2, 3, 5, 10 and 20 layers, shows that the topology of the phase 
diagrams is not affected by the increasing film thickness, for a fixed exponent a value. Hence, 
the numerical results established in this work are done for a film with a thickness = 10 
layers. 

In order to outline the effect of the parameter a at a fixed film thickness, we plot in Figs. 2a, 
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2b, 2c and 2d the corresponding phase diagrams for several a values: 0.0, 0.5, 1.0 and 2.0. 
In the case a — 0.0, Fig. 2a, a constant crystal field is applied on each layer and the corre- 
sponding transition is found at A^/J = 3 — = 2.9 for very low temperatures. 
Increasing a. values, it is seen that large values of A^/ J are needed to disorder the film. This 
result is summarised in Figs. 2b, 2c and 2d. Indeed, due to the effect of the equation law (4), 
separate layering transitions are found. For non null but very low temperatures, the surface 
crystal field values found in the phase diagrams are exactly those predicted by the ground state 
study: A^/J = 3(/c)" for a layer order k with k < k^. For example, in Fig. 2c for a — 1.0, 
the first transition is found at A^/J = 3.0, the second transition at A^/J = 6.0, the third 
transition at A^/ J = 9.0 and so on. The other numerical values of the surface crystal field 
corresponding to a = 0.5 in Fig. 2b, and a = 2.0 in Fig. 2d are exactly those predicted by the 
ground state phase diagram. Furthermore, Figs. 2b, 2c and 2d show that the remaining N — ko 
transitions, not seen for T = 0, arise for higher temperature values with separate tricritical 
points. More over, these results ensure the existence of a critical value a^k) above which a 
layer k disorders at a very low temperature. The reentrant phenomena, seen in these figures 
for a given layer k, is caused by the competition of the temperature and an inductor magnetic 
field created by the deeper layers. Indeed, when the thermal fiuctuations become sufficiently 
important, the magnetisation of some spins, of deeper layers, becomes non null (+1 or -1). This 
leads to the appearance of an inductor magnetic field. This magnetic field is responsible of the 
ordered phase seen for the layer k. This argument can also explain the absence of the reentrant 
phenomena for the last layer (s), once the magnetisation of the remaining layer(s) — A; is not 
sufficient to create an inductor magnetic field. It is worth to note that the reentrant phenom- 
ena is always present for the layers k, {k < ko), and the corresponding tricritical points Cj are 
located at a constant temperature. An interesting question arises at this stage: what is the 
effect of the parameter a on the surface of the film (fc = 1)? To clarify this point, we illustrate 
in Fig. 3a the profiles of the critical temperature for several values of o; as a function of As/ J. 
It is found that increasing a values leads to important reentrants phenomena. Indeed, for large 
values of a, the layer with k > 1 needs higher values of the surface crystal field to disorder. 
It is also found that, except the special case a — 0.0, the tricritical point is not affected by 
variations in a values. This surface tricritical point is located at (A^/J = 3.0, Tc/ J = 0.78) 
for a 7^ and {As/ J = 2.9, TJ J = 1.9) for a = for a film with a thickness = 10 layers. 
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The effect of increasing a on a medium layer k = N/2 = 5 is illustrated in Fig. 3b. The same 
topology is found, but higher values of A^/J are needed to overcome the disorder of this layer. 
On the other hand, for o; — > oo the studied system is equivalent to a configuration with 
As/J{k = 1) = cste and As/J{k > 1) = 0.0: the applied crystal field is present only at the 
surface and absent on the bulk. The corresponding phase diagram is illustrated by Fig. 4. 
Only two transitions are present in this case: 0^° ^ DO^, occurring for low temperatures at 
Ag/J — 3.0, and DO^ D^^ occurring at a constant temperature T/J = 3.899 and any non 
null surface crystal field value. In order to complete this study, we illustrate in Fig. 5, the ther- 
mal behaviour of the surface magnetisation for A^/ J = 3.2 and selected values of a: 0.0, 0.2, 
0.5, 1.0 and 2.0. At this surface crystal field value, the increasing a values effect is to increase 
both the surface magnetisation amplitude and the critical temperature needed to disorder this 
layer. This is a consequence of the factor A;" effect for A; > 1, as it has been anticipated in the 
discussion above. 



4 Conclusion 

We have studied the effect of a variable crystal field according to the law Eq. (4), and found 
order-disorder layering transitions of the Blume-Capel Ising model using mean field theory. We 
found a reentrant phenomena for each layer k {k < ko, see text), and a tricritical point at a fixed 
temperature with a vertical first order line. We established the temperature-crystal field phase 
diagrams and found that the last layer tricritical point depend strongly on the small values of 
the exponent a for a fixed film thickness. The thermal behaviour of the surface magnetisation 
has also been investigated. 
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Figure Captions 



Figure 1.: The ground state phase diagram in the plane (J, A^) for a film thickness N. The 
notations D^O^~^ are defined in the body text. The calculations presented in this manuscript 
correspond to a film thickness N — layers. 

Figure 2.: The critical temperature behaviour as a function of the surface crystal field for 
several values of the exponent a : a — 0.0, h) a — 0.5, c) a — 1.0 and d) a — 2.0. 
The first-order transition hues (dashed hue) are connected to the second-order transition lines 
(continuous line) by a tricritical point. 

Figure 3.: The critical temperature behaviour of the first layer k = 1 and the medium layer 
A; = 5, as a function of the surface crystal field A^ for selected values of the parameter a: 
0.0, 0.2, 0.51.0 and 2.0. a) First layer A; = 1. b) Medium layer k = N/2 = 5. 

Figure 4.: The dependence of the critical temperature as a function of the surface crystal field 
A5 for a — 00. Only two transitions are present in this case: 0^° ^ DO^ and DO^ ^ D^^. 

Figure 5.: The thermal behaviour of the surface (layer A; = 1) magnetisation for selected values 
of a: 0.0,0.2,0.5,1.0 and 2.0. 
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